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' _h' ' Abstract 

Mh' The aim of this article is to propose a mathematical model describing the elec- 

ta , tronic structure of crystals with local defects in the framework of the Thomas-Fermi- 

von Weizsacker (TFW) theory. The approach follows the same lines as that used in 
E. Cances, A. Deleurence and M. Lewin, Commun. Math. Phys., 281 (2008), pp. 
129-177 for the reduced Hartree-Fock model, and is based on thermodynamic limit 
arguments. We prove in particular that it is not possible to model charged defects 
within the TFW theory of crystals. We finally derive some additional properties of the 
^ , TFW ground state electronic density of a crystal with a local defect, in the special case 

\mml • when the host crystal is modelled by a homogeneous medium. 

• ' 1 Introduction 

O, 

("^ I The modelling and simulation of the electronic structure of crystals is a prominent topic in 

solid-state physics, materials science and nano-electronics [13 1 H^ j [T9]. Besides its impor- 
tance for the applications, it is an interesting playground for mathematicians for it gives 

k>( , rise to many interesting mathematical and numerical questions. 

^ ' There are two reasons why the modelling and simulation of the electronic structure of 

crystals is a difficult task. First, the number of particles in a crystal is infinite, and second, 
the Coulomb interaction is long-range. Of course, a real crystal contains a finite number of 
electrons and nuclei, but in order to understand and compute the macroscopic properties 
of a crystal from first principles, it is in fact easier, or at least not more complicated, to 
consider that we are dealing with an infinite system. 

The first mathematical studies of the electronic structure of crystals were concerned 
with the so-called thermodynamic limit problem for perfect crystals. This problem can 
be stated as follows. Starting from a given electronic structure model for finite molecular 



systems, find out an electronic structure model for perfect crystals, such that when a cluster 
grows and "converges" (in some sense, see [4J) to some 7^-periodic perfect crystal, the ground 
state electronic density of the cluster converges to the 7^-periodic ground state electronic 
density of the perfect crystal. 

For Thomas-Fermi like (orbital-free) models, it is not difficult to guess what should be 
the corresponding models for perfect crystals. On the other hand, solving the thermody- 
namic limit problem, that is proving the convergence property discussed above, is a much 
more difficult task. This program was carried out for the Thomas-Fermi (TF) model in |l6j 
and for the Thomas- Fermi- von Weizsacker (TFW) model in 0]. Note that these two models 
are strictly convex in the density, and that the uniqueness of the ground state density is 
an essential ingredient of the proof. The thermodynamic limit problem for perfect crystals 
remains open for the Thomas-Fermi-Dirac-von Weisacker model, and more generally for 
nonconvex orbital- free models. 

The case of Hartree-Fock and Kolm-Sham like models is more difficult. In these models, 
the electronic state is described in terms of electronic density matrices. For a finite system, 
the ground state density matrix is a non-negative trace-class self-adjoint operator, with trace 
A^, the number of electrons in the system. For infinite systems, the ground state density 
matrix is no longer trace-class, which significantly complicates the mathematical arguments. 
Yet, perfect crystals being periodic, it is possible to make use of Bloch-Floquet theory and 
guess the structure of the periodic Hartree-Fock and Kohn-Sham models. These models 
are widely used in solid-state physics and materials science. Here also, the thermodynamic 
limit problem seems out of reach with state-of-the-art mathematical tools, except in the 
special case of the restricted Hartree-Fock (rHF) model, also called the Hartree model in 
the physics literature. Thoroughly using the strict convexity of the rHF energy functional 
with respect to the electronic density, Catto, Le Bris and Lions were able to solve the 
thermodynamic limit problem for the rHF model [5]. 

Very little is known about the modelling of perfect crystals within the framework of the 
A^-body Schrodinger model. To the best of our knowledge, the only available results |8l |12] 
state that the energy per unit volume is well defined in the thermodynamic limit. So far, 
the Schrodinger model for periodic crystals is still an unknown mathematical object. 

The mathematical analysis of the electronic structure of crystals with defects has been 
initiated in [1] for the rHF model. This work is based on a formally simple idea, whose 
rigorous implementation however requires some effort. This idea is very similar to that used 
in [6l [ini [11] to properly define a no-photon quantum electrodynamical (QED) model for 
atoms and molecules. Loosely speaking, it consists in considering the defect (the atom or 
the molecule in QED) as a quasiparticle embedded in a well-characterized background (a 
perfect crystal in our case, the polarized vacuum in QED), and to build a variational model 
allowing to compute the ground state of the quasiparticle. 

In [Ij, such a variational model is obtained by passing to the thermodynamic limit in 
the difference between the ground state density matrices obtained respectively with and 



without the defect. In order to avoid additional technical difficulties, the thermodynamic 
limit argument in [Ij is not carried out on clusters (as in jH [5]), but on supercells of 
increasing sizes. Recall that the supercell model is the current state-of-the-art method to 
compute the electronic structure of a crystal with a local defect. In this approach, the defect 
and as many atoms of the host crystal as the available computer resources can accomodate, 
are put in a large, usually cubic, box, called the supercell, and Born-von-Karman periodic 
boundary conditions are imposed to the single particle orbitals (and consequently to the 
electronic density). The limitations of the supercell methods are well-known: first, it gives 
rise to spurious interactions between the defect and its periodic images, and second, it 
requires that the total charge contained in the supercell is neutral (otherwise, the energy 
per unit volume would be infinite). In the case of charged defects, the extra amount of 
charge must be compensated in one way or another, for instance by adding to the total 
physical charge distribution of the system a uniformly charged background (called a jellium). 
It is well-known that this procedure generates unphysical screening effects. Other charge 
compensation methods have been proposed, but none of them is completely satisfactory. 
Note that the above mentioned sources of error vanish in the thermodynamic limit, when 
the size of the supercell goes to infinity: both the interaction between a defect and its 
periodic images and the density of the jellium go to zero in the thermodynamic limit. 

The variational model for the defect, considered as a quasiparticle, obtained in [T] has a 
quite unusual mathematical structure. The rHF ground state density matrix of the crystal 
in the presence of the defect can be written as 

7 = 7per + Q 

where 7pgj, is the density matrix of the host perfect crystal (an orthogonal projector on 
L^(]R^) with infinite rank which commutes with the translations of the lattice) and Q a self- 
adjoint Hilbert-Schmidt operator on L^(M^). Although Q is not trace-class in general [2], 
it is possible to give a sense to its generalized trace 

Tro(Q) := Tr(Q++)+Tr(Q— ) where Q++ := {l-^l^^)Q{l-^l^^) and Q— := 7°erQ7per 

(as 7ppj, is an orthogonal projector, Tr = Trg on the space of the trace-class operators on 
L^(M )), as well as to its density pq. The latter is defined in a weak sense 



yW £ C^{R^), TroiQW) = I pqW. 



The function pq is not in L"'^(]R'^) in general, but only in L^(M'^)nC, where C is the Coulomb 
space defined by (J12p . An important consequence of these results is that 



in general, the electronic charge of the defect can be defined neither as Tr(Q) nor as 

/r3P 



• it may happen that pq G L"'^(M'^) but Ttq{Q) ^ f^a PQ (while we would have pq E 
L-^(M^) and Tro(Q) = Tr((3) = La Pg if Q were a trace-class operator). In this 
case, Tiq{Q) and Jj^s pq can be interpreted respectively as the bare and renormalized 
electronic charges of the defect [2] . 

The reason why, in general, Q is not trace-class and pq is not an integrable function, is 
a consequence of both the infinite number of particles and the long-range of the Coulomb 
interaction. 

Note that, still in the rHF setting, the dynamical version of this variational model is 
nothing but the random phase approximation (RPA), widely used in solid-state physics. 
The well-posedness of the nonlinear RPA dynamics, as well as of each term of the Dyson 
expansion with respect to the external potential, is proved in [3]. 

As far as we know, the mathematical study of the electronic structure of crystals with 
local defects has not been completed for the Thomas- Fermi- von Weizsacker model |14) . This 
is the purpose of the present work. The article is organized as follows. In Section [21 we 
present the periodic TFW model used in condensed phase calculations. After recalling the 
mathematical structure of the TFW model for perfect crystals (Section 13. ip . we propose 
a variational TFW model for crystals with local defects (Section 13. 2p . We prove that this 
model is well-posed and that the nuclear charge of the defect is fully screened, in a sense 
that will be precisely defined. In Section 13.31 we provide a mathematical justification of 
the model introduced in Section [3] based on bulk limit arguments. In Section [3.41 we focus 
on the special case when the host crystal is a homogeneous medium, that is when both the 
nuclear and electronic densities of the host crystal are uniform (and opposite one another 
to prevent Coulomb blow-up) . The technical parts of the proofs are gathered in Section [H 

Note that the screening effect has already been studied in the context of the Thomas- 
Fermi model in |16j . in the case when the host crystal is a homogeneous medium. 



2 The periodic Thomas- Fermi- von Weiszacker model 

In this section, we describe the Thomas-Fermi-von Weiszacker (TFW) model with Born- 
von Karman periodic boundary conditions, used to perform calculations in the condensed 
phase. In Section 13.31 we will use this periodic model to pass to the thermodynamic limit 
and construct a rigorously founded TFW model for crystals with local defects. 

Let 7^ be a periodic lattice of R^, 7^* the associated reciprocal lattice, and P the simu- 
lation cell. If for instance 7?. = aZ^ (cubic lattice of size a), then 7^* = -^^'^ and possible 
choices for P are P = (0, a]^ or P = (— §, f]^- Let also P* be the first Brillouin zone of the 
lattice TZ (or in other words, the Wigner-Seitz cell of the reciprocal lattice TZ*). 

We introduce the usual 7^-periodic L^ spaces defined by 

LP,,(P) := {v G Ll^iM.^) \ v 7l-periodic} , 



and endow them with the norms 



V rP 



LlcriT) •= ( / 1^1^ ) for l<p<oo and ||i'||Lgg^(r) := ess-sup|t;| 

In particular, 

ll^llL2^,(r) = (v, v)Ll^(r) wliere {v, w)^^^^^^) := j^ vw. 
Any function v £ Lpgj.(r) can be expanded in Fourier modes as 

„ik-x 



sr-^ e 1 f 

(x) = 2_] ck{v)—— where Ck{v) = — — / ^ 

1,^'D* II II •'T 



v(x)e "-^'"^ dx. 
ken* 



The convergence of the above series holds in Lp^j.(T, C), the space of locally square integrable 
7?.-periodic C-valued functions. 

For each s G M, the "/^.-periodic Sobolev space of index s is defined as 

v{x) = Y, Cfc W^ I E (1 + l^l')1cfc Wl' < oo, VA: G 7^^ c_fc = c^ , 
fceT?.* ' ' ken* ) 

and endowed with the inner product 

The condition VA; G 7^*, c__fc = c^ implies that the functions of Hp^^(T) are real- valued. 
Recall that H^^^{T) = Ll^.iT), {■,-)H0^^{r) = {■,-)Ll^^{r), 

^per(r) = {ve Ll,,{r) I V^ G (L2,,(r))'} , {v, w)Hi^^^r) = Jvw + JvvVw, 



and (Fp-?(r))' = ^-,r(r). 

We also introduce the 7^-periodic Coulomb kernel G-n defined as the unique function of 

^pcrv 



Lpgj,(r) solution of the elliptic problem 



-AG7^ = 47rK]5fc-|rrM 
\ken ) 

Gti 7^-periodic, minGT^ = 0. 

It is easy to check that 

Ak-x 






ken*\{o} 



The 7^-periodic Coulomb energy is then defined for ah / and g in Lpj,j,(r) by 
DnifiO) = / I Gn(x-y)f{x)g{y)dxdy 



r JT 



Gn]co{f)co{g)+ V Tr72^k{f)ck{g) 
■>^ ^ fcg7^•\{o} ' ' 

{Gn ^n f){y)9{y) dy= {Gn -kn g){x)f{x) dx, 
r Jr 

where -k-ji denotes the 7^-periodic convolution product: 

v(/,5) G ^per(r) X Ll,,{T), (/ *n g){x) = I f(^- y)9{y) dy = J f{y)g{x - y) dy. 

Let p™"^ be a function of Lpgj.(r) modelling a 7^-periodic nuclear charge distribution (or 
the effective charge distribution of a pseudopotential describing a 7^-periodic distribution of 
nuclei and core electrons) . The corresponding 7^-periodic TFW energy functional is defined 
on i7pg^(r) and reads 

^TFW(pnuc^^) =Cyv [ \Vv\^ + Ctf / 1^1^°/^ + ^Dnip''''^ " ^'^ /o'^"' - v^), (1) 
Jr Jr 2 

where 

Ctf = — (37r ) ' (Thomas-Fermi constant) and Cw > 
o 

(several values for Cw have been proposed in the literature, see e.g. [7]). From a physical 

viewpoint, p = v^ represents the electronic density (or the electronic density of the valence 

electrons if the core electrons are already incorporated into p^^^). The first two terms of 

E'^^ {p™^'^ ,v) model the kinetic energy per simulation cell and the third term the Coulomb 

energy of the total 7^-periodic charge distribution /9*°* = p™^^ — v^ . 

The electronic ground state with Q electrons in the simulation cell is obtained by solving 

the minimization problem 

/7^(p""^Q) = inf |i5;™(^nuc^^)^ ^ ^ ^i^^(r), J/ = q]- (2) 

For the sake of simplicity, we assume that the nuclear charge density is in Lpg^(r). This 
allows us to gather all the Coulomb interactions in a single, non-negative term (the third 
term in the right hand side of ([T])). On the other hand, this excludes point-like charges 
represented by Dirac measures. As often in this field, it is however easy to extend our 
analysis to point-like nuclei, by splitting the Dirac measure 5q as 5q = {5q — (p) + (j) where 
(/> is a radial function of C^(M^) such that Jj^g </> = 1 and Supp(0) small enough. 

The following result is classical. We will however provide a proof of it in Section |4] for 
the sake of completeness. 



Proposition 2.1. Let p^"^"" G -^pcr(r) and Q > 0. 

1. Problem ([^ has a minimizeru such that u G Hp^j.{r) "^ C^(M^)nL°°(M^) andu > 
in M^. The function u satisfies the Euler equation 

- Cw^u + \ctW^ + {Gn *7^ («' - P°"')) n = e^u, (3) 

where ep «s i/te Lagrange multiplier of the constraint J-^u^ = Q. 

2. Problem ^ has exactly two minimizers: u and —u. 

As a consequence of Proposition 12.11 the ground state electronic density is always 
uniquely defined in the framework of the periodic TFW model. 

3 The Thomas-Fermi- von Weiszacker model for crystals 

We now focus on the special case of crystals. More precisely, we consider two kind of 

systems: 

• a reference T^i-periodic perfect crystal with nuclear distribution 

Pper ^ -^per(ri), 

where Fi is a unit cell for IZi] 

• a perturbation of the previous system characterized by the nuclear distribution 

C denoting the Coulomb space defined by (J12p . 

3.1 Reference perfect crystal 

It is shown in [4j that the ground state electronic density /Opgj. of a crystal with nuclear 
charge distribution /Op^Jr G L'l^^^iTi) can be identified by a thermodynamic limit argument. 
It is given by p'i^^ = jUperP where Upg^ > is obtained by solving the minimization problem 

/7^,(p-^Z)=inf|i?TFW(^nuc^^)^ ^^^l^^(ri), ^ V^ = Z^, (5) 

where 

Z= f Pll'i. (6) 



Note that problem ([5]) has a unique solution (up to the sign) for any value of Z. The correct 
value of Z given by ([6]) is obtained in [4J by a thermodynamic limit argument. As expected, 
this value implies the charge neutrality condition 

{p7^ - pI..) = 0. (7) 



The unique non- negative solution Up^j. to ([S])-® satisfies the Euler equation 

- CwA-U^pj. + -CTF(Pper)^ ^^per + (^-^^i *-?^i (Pper " Ppe?)) ""per = ^P'^per' (8) 

where e^, the Lagrange multiplier of the charge constraint, which is uniquely defined, 
is called the Fermi level of the crystal. From ([7|, we infer that the Coulomb potential 
V?gj. = G-ji^ i^Ti^ (Pper ~ Pper ) ^^ ^^ Unique solution in }l\^^^{V\) to the T^i-periodic Poisson 
problem 

•^per i-p I \rpcr rpcr/ ' 



Fp^er 7^l-periodic, / V^^, = 0. 



By elliptic regularity, V^^,. G H^crO^i) ^ C'^{R^) D L'^{M^). Using Proposition EH we 
obtain that u^^^ G C^{R^) D L°°Im^), and that -u^g^ > in M^. We thus have the following 
bounds, that will be useful in our analysis: 

30 < m < M < +00 s.t. Vx G M^, m< ul^j.{x) < M. (9) 



Let us denote by H^^^ the periodic Schrodinger operator on L (M ) with domain H {. 
and form domain //^(M^) defined by 



^V G F2(M3), Hl,,V = -CwAv + -CTFiplerf/^V + ^ 



V. 



It is classical (see e.g. [TSj) that H^^^j. is self-adjoint and bounded from below, and that its 
spectrum is purely absolutely continuous and made of a union of bands. For convenience, 
we will use the abuse of notation consisting in denoting by H9^^^v the distribution 



per'- 

5 



fpVz; := -C^Av + -CMplrf^'v + V^ 



which is well-defined for any v G Lj'q^(M^), and belongs to H ^(R^) if t; G H^(R^) and to 
H~Jj.{T) if f G //pgj,(r). We can thus rewrite equation ([8]) under the form 

-"per "per " ^F^pcr- \^^} 

Using the fact that u^^j. > in M^, it is easy to see that e^ is in fact the minimum of the 
spectrum of the periodic Schrodinger operator i:^per (that is the bottom of the lowest energy 
band). As a consequence, 

V^ G H\r'), {{H^,, - e^)^,f)H-i(R3),,^i(M3) > 0. (11) 



3.2 Crystals with local defects 

We now consider a crystal with a local defect whose nuclear charge distribution is given by 
(J4]). It is convenient to describe the TFW electronic state of this system by a function v 
related to the electronic density p by the relation 

We denote by C the Coulomb space defined as 

C := {/ G 5'(R3) I / G LlM^), I • ri/(.) G L\R')] , (12) 

where / is the Fourier transform of /, normalized in such a way that ||/||l2(r3) = ||/||l2 
for all / G L'^{M.^). Endowed with the inner product 



D{f,g) -.= 47^ [ 



f{k)g{k) 



dk, 



C is a Hilbert space. It holds L^/^{R^) C C and 



°"^ - u-^l 



\X — X 

Denoting by 

JiRS JlR3 2 

the TFW energy functional of a finite molecular system in vacuo with nuclear charge p^^''^ , 
we can formally define the relative energy (with respect to the perfect crystal) of the system 
with nuclear charge density Pp^Jr + u and electronic density p = {u^^^ + v)'^ as 

E^^'^ipll^ + ., ul, + v)- E^^^ipl^^, u%,) 

+\d {2ul,,v + v^-v, 2ul,,v + v^-v)- f vV^,, + eU, (13) 

where 

.0 , „,|2 I ,0 |2 



per 



+ v\^-\ul,f). (14) 



Of course, the left-hand side of (|13p is a formal expression since it is the difference of two 
quantities taking the value plus infinity. On the other hand, the right-hand side of (J13p is 
mathematically well-defined as soon as g is a fixed real number and v G Q+, where 

Q+ := {v e HHM.^) I V > -u%,, u%,v G C} . 



The set Q+ is a closed convex subset of the Hilbert space 

endowed with the inner product defined by 

{v,w)q := {v,w)hH^3) +L»(UpcrW,n°cr^«)- 



This formal analysis leads us to propose the following model, which will be justified 
in the following section by means of thermodynamic limit arguments: the ground state 
electronic density of the perturbed crystal characterized by the nuclear charge density Q 
is given by 

Pu = (Uper + Vuf, 

where v^ is a minimizer of 

r = mf{£''{v),veQ+} (15) 

with 

5 
3' 

+-D {2ul^^v + v'^ -u, 2ul^^v + v'^ -u). (16) 

The following result, whose proof is postponed until Section HI shows that our model is 
well-posed. 

Theorem 3.1. LetuGC. Then, 

1. Problem ( Ii5j) has a unique minimizer v^, and there exists a positive constant Cq > 
such that 

^y^C, lb.||Q<Co(||i^||c + lkl|^). (17) 

The function v^, satisfies the Euler equation 



+Ctf I ( \u%, + v\''/' - \u%,r/'' - ^-\ulJ'/H2u%,v + v') 



tuO ^0 N,,, I ^r* (U,0 _i_„ |7/3 u,0 |7/3 u,0 |4/3„ 



5 

per ^iij^v I ■;j'-'ir \ | ""per ' "v\ I ""perl I ""perl 

+ [{2ul,,v, + vl-v)^\-\'^) {ul,, + V,) = 0. (18) 

2. Let us denote by p^ = v — {2u^f,^Vy + v^) the total density of charge of the defect and 
by ^^ = p^ -k\ ■ \^^ the Coulomb potential generated by p^. It holds Vu G H'^iM.^), 
$0 G L2(M3) and 

hm-^ / \p^{k)\dk = 0. (19) 

r^O \Br\ J Br 

10 



3. Any minimizing sequence {vn)neN for il5]) converges to Vy weakly in H^(^) and 
strongly in Lj^^(M^) for 1 < p < 6. Besides, (Wper'^n,)neN converges to v^^^Vy 
weakly in C. 

For any q G M, there exists a minimizing sequence (fn)nGN for U5\) consisting of 
functions of Q^ n L^(R^) such that 

VnGN, / {\ul,,-,VrS^-\ul,,\^)=q. (20) 



We conclude this present section with some physical considerations regarding the charge 
of the defect. 

Remark 3.1. Let v E L1(M3) n L'^i^). Assuming that v^ E L^{R^) D L^{R^) (a property 
satisfied at least in the special case of a homogeneous host crystal, see Section \'J.4\), then 



p1 E C''(M^) and I119\) simply means that the continuous function fp^ vanishes at k = 0, or 
equivalently that 

' pI = 0. (21) 



The property il9\) means that is a Lebesgue point of p^ and that the Lehesgue value of 
p^ at is equal to zero. It can therefore he interpreted as a weak form of the neutrality 
condition \21\) . also valid when p^ ^ L-'^(M^). The fourth statement of Theorem \3.1\ implies 
that there is no way to model a charge defect within the TFW theory: loosely speaking, if 
we try to put too many (or not enough) electrons in the system, the electronic density will 
relax to (Wper + ^i/)^ '^'^^ the remaining (or missing) q — J^^ v electrons will escape to (or 
come from.) infinity with an energy ep. 

3.3 Thermodynamic limit 

The purpose of this section is to provide a mathematical justification of the model (|15p . 
Consider a crystal with a local defect characterized by the nuclear charge distribution 

pnuc^^nuc^j^ ^-^j^ z^EL^(M3)nL2(M3). (22) 

In numerical simulations, the TFW ground state electronic density of such a system is 
usually computed with the supercell method. For a given L E N large enough, the supercell 
model of size L is the periodic TFW model ([2]) with 

n = nL:=Lni, r = ri:=Lri, p'^"= = p^»^= + i.i, Q = zL'^ + q, (23) 

where 

M^) = Yl (xri,i^)(x-z), 

11 



XFl : ]R^ — ^ IK denoting the characteristic function of the simulation cell F^. Note that v^ 
is the unique T^j^-periodic function such that i^l\tl — ^Ir^,- In practice, L is chosen as large 
as possible (given the computational means available) to limit the error originating from 
the artificial Born-von Karman periodic boundary conditions. 



It is important to note that Up^j, is the unique minimizer (up to the sign) of the supercell 
model of size L for p^'^^ = p^'^ and Q = ZL^, whatever L G N*. Reasoning as in the 
previous section, we introduce the energy functional 

+C7TF l^ {\nlr + vl\''/' - Wir/' - lKj'/'{2n%.v, + .i)) 

+ 2^^L (2'Uper^L + vl- lyi, 2ulej.VL + vl- vl) , (24) 



which is such that 



Jtt. 



with 



/ (Kr + VL?- K/) = / {2U%VL + vl). (26) 



While (|13p and (J14p are formal expressions, (j25j) and (j26p are well-defined mathematical 
expressions. The ground state electronic density of the supercell model for the data defined 
by (p3|) is therefore obtained as 

Pl = (""per + Vu,q,L) 

where v^^q^i is a minimizer of 

11" = inf S^S'iivL), VL G Q+,L, ^ i2u'^erVL + ^1) = oj , (27) 

Q+^L denoting the convex set 

Q+,L = {VL(^ -f^per(rL) | Vl > -U^^r] ■ 

We also introduce the minimization problem 

rL = -^nf{£UvL),VLeQ+,L}, (28) 

in which we do not impose a priori the electronic charge in the supercell. 



12 



Theorem 3.2. Let v £ L^{R^) n L 

1. Thermodynamic hmit with charge constraint. For each q £ M and each L G N*, 
the minimization problem i27\) has a unique minimizer Vu^q^L. For each q £ M, the 
sequence (Wi/,g,L)LgN* converges, weakly in H^^^{M.^), and strongly in L^^^{W^) for all 
1 < p < Q, towards v^, the unique solution to problem ( [J5|) . For each g E M and each 
L S N*, Vv^q^L satisfies the Euler equation 



(ttO ^0 \^, I '-'r' lla.O _L ,, |7/3 U,0 |7/3 u,0 |4/3 



+ ((2Uper^^i.,g,L + ^^^,g,L " T^l) -^TZl G-Rl) (^por + Vv,q,L) = /^z.,g,L (Wper + Vv,q,L), (29) 

where fJ-u,q,L £R is the Lagrange multiplier of the constraint Jp {'^Up^j.Viy^q^L-'rv'^ i) = 
q, and it holds hm fiu,q,L = for each g G M. 

2. Thermodynamic limit without charge constraint. For each L G N*, the minimization 
problem \28\) has a unique minimizer v^^l. It holds 

fuO ^0\^, i^n /'lo.O I „ |7/3 |,,0 |7/3 |,,0 |4/3,, 

l^per - ^F)Vu,L + -C-TF (^Pper + Vu,l\ " Pperl " Pperl ^^^.i 

+ ((2<r^^,L + vl,L - ^l) ^TIl Gul) «r + Vu,l) = 0. (30) 

The sequence {vu,L)LeN* o-lso converges to Vy, weakly in Hyoc^^)> ^^'^ strongly in 
Lf^jR^) for all 1 < p < 6. Besides, 






3.4 The special case of homogeneous host crystals 

In this section, we address the special case when the host crystal is a homogeneous medium 
completely characterized by the positive real number a such that 

Vx G M^, pl'^^ix) = pI^^{x) = a^ and ^^^^(x) = a. (31) 

In this case, analytical expressions for the linear response can be derived, leading to the 
following result. 

Theorem 3.3. Assume that \'J1\) holds. For each v £ C, the unique solution Vy to Iil5\) can 
be expanded as 

Vu = g-ku + ¥2(2^) (32) 

where g £ L^(M'^) is characterized by its Fourier transform 

1 Ana 



9{k) 



(27r)3/2 Cw\k\^ + f CTFa4/3|fc|2 + 8to2 ' 



9 
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and where r2(y) G L^{^^). For each u G L^{E?) n C, it holds v^ G L^(m?) n L'^{M?) and 

f {u - i2u%,v^ + vD) = 0. 

The first term in the right hand side of ()32p is in fact the linear component of the 
application ly i-^ Vi,. The second term gathers the higher order contributions. 



Proof. In the special case under consideration, the Euler equation (jlSp also reads 

20 

-Cw^Vu + —CTFa'^^^v„ + 2a^ (w,, • | • |~^) = a (i/* | • |"^) - a (w^ * | • |"^) + k^, (33) 
9 



where 



^Ctf (\a + f,r/=^ - a^/=^ - ^a^/^^'.) + ((^^ - '^av, -vl)ir\. \-') v,. 



3 

We therefore obtain (|32p with 

r2i'^) = -g-kvl + h-kKiy, 
the convolution kernel h being defined through its Fourier transform as 

1 \k\^ 



h{k) 



(27r)3/2 Cw|A;|4 + fCTYa^/^\k\^ + Svra^ ' 



It follows from the second statement of Theorem 13.11 and Lemma 14.11 below that k^, G 
L"^(M3). The proof will therefore be complete as soon as we have proven that g G L"'^(M3) 
and h G ^^(M^). In fact, we will prove that any even-tempered distribution / G 5'(M'^) 
whose Fourier transform is a function of the form 



|fc|r(H,-|) 



where q and r are polynomials of the real variable satisfying deg(g) < deg(r), r > in M+, 
g(0) = 0, q"{Q) = and r'(0) = (which is the case for both g and h), is in L^(M'^). Indeed, 
/ G L2(R3) and a simple calculation shows that 



/(x) = a/^^ r"(^)(t)sin(|x|t)di 
V TT \x\ Jq \rJ 
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Therefore, there exists C G M_(_ such that 

C 



1/(^)1 < 



I 19 I I I ^ 

almost everywhere in M^, which proves that / G L^(M^). D 

Remark 3.2. For a generic z/ G C, the function Vy, hence the density 2aVy + v^, are not 
in L^(M^). This follows from the fact that the nonlinear contribution r2{y) is always in 
-L^(M^), while the linear contribution g -k u is not necessarily in L^{M.^) since its Fourier 
transform, 



{g-kv){k) = ^ ,,|4 20r^ ./.I, I. . . — ^^W 



is not necessarily in L^ 

4 Proofs 

This section is devoted to the proofs of Proposition 12.11 Theorem 13.11 and Theorem 13.21 
In the sequel, we set 

Ctf = 1 and Cw = 1 (in order to simplify the notation). 

4.1 Preliminary results 

We first state and prove a few useful lemmas. Some of these results are simple, or well- 
known, but we nevertheless prove them here for the sake of self-containment. 

Lemma 4.1. For all < m < AI < oo and all 7 > 2, there exists C G M+ such that for 
all m < a < M and all b > —a, 

(7 - l)a^-2^2 < (^ ^ 5)7 _ ^7 _ ja^-^b <C{1 + 16^-2) 52. (34) 

Proof. Let (j){t) = (a + tb^ . It holds for all t G (0, 1), (/>'(t) = 7(0 + tb^-'^b and (/)"(t) = 
7(7 — l){a + tb)'^~'^b'^. Using the identity 

(^(1) - <f>{0) - (/.'(O) = / (1 - t)(t>"{t) dt, 
Jo 

we get 

(a + by -a^ - ja^-% = 7(7 - l)^^ f (1 _ t)(a + t6)7-2 ^^^ 

Jo 

We obtain dMI using the fact that for all t £ [0,1], a{l - t) < a + tb < M + \b\. D 
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Lemma 4.2. Let u e C and v e Q+ n H'^{R^) such that v > -Uper in M.^. For all e > 
and g G M, there exists v^ G Q+ H C^(M^) such that 

I (2u° Ve + v^,) = q and {S^iv,) - £''{v)\ < e. 

JR3 



Proof. Let e > 0. As the functions of H'^(U.^) are continuous and decay to zero at infinity, 
there exists 6 > such that 

Vx G M^ v{x) > -ul^,{x) + 5. (35) 

For all R > 0, let Bji be the ball of M'^ centered at zero and of radius R. For 77 > 0, we 
define 

where T is the Fourier transform and J^^^ the inverse Fourier transform. Clearly, v^ G 
H (R.^) ^-7- C^(M^) and Up^j-v"^ G C. Li addition, when r] goes to zero, {v'^)n>o converges to 
V in H'^{R^), hence in L°°(M'^), and {Uper''j'^)'q>o converges to Up^j.v in C. The function S'^ 
being continuous on Q, this implies that there exists some r^o > such that 

t;* G Q+ n C\R^) and IS'^iv'^^) - S^v)] < e/4. 

Let X be a function of C^(M'^) supported in B2, such that < x(") ^ 1 ^iid X = 1 iii 
i?i. For n G N*, we denote by Xn{') = x('^~^") ^-nd by v^°'^ = XnV^°- For each n G N*, 
yVo,n g Q_|_ p (7^(1^3^ g^jjj ^]^g sequence (z;''°'"')„gN* converges to v'^° in Q when n goes to 
infinity. Hence, we can find some uq > such that 

yrio,no e g^ n C2(m3) and |£:^(t;*'"o) - £:'^(t;'?o)| < e/4. 
Let 

90 



/ (2tx° ^""'"^ + (^''°'"°)^) and q^ = q-qo. 

Jr3 



If (71 = 0, f "^ = t;'?0'"o fulfills the conditions of Lemma 14.21 Otherwise, we introduce for m 
large enough the function Vm defined as Vm = tmXmUper where tm is the larger of the two 
real numbers such that 

/ (,'^^per'^m ' '^m) ~ 2lm / XraPper ~^ '^m I XmPpcr ~ 9l- 

Jm.3 "^ Jk.3 "^ Jm.3 

A simple calculation shows that tm ~ — 5i|ri|Z^ ( / x] m , and that 

m^oo 2 \JjK3 y 

lim £:0(z;„) = 0, 
16 



so that there exists tuq £ N* such that Vm G Q+ n Cc(M^) and < £^{vmo) < e/4. Let 
us finally choose some Ri G TZi \ {0} and introduce the sequence of functions (wmoj? )peN 
defined by 

<o?P°(-)=^''""°(-) + ^™o(--P^l)- 

For p large enough, Vmoji belongs to Q+ n C^(M^) and satisfies 

/ (2<e.<rp° + «7f) = 9- 

Besides, 

< e/4 + \D{2u%,v^^'^'> + (^^"'"0)2 - j., {2u%,Vmo + <)(• " P^i))| • 

As 

^l|m Z)(2«0,,^^0'-» + (^^"'"0)2 - z., (2^0,,z;„„ + vl^){- - pR,)) = 0, 

there exists some pQ G'N such that 

\D{2u%,v^'^^^^ + (v^^'-'^f - V, {2u%,vm, + vlj{- - pRi))\ < 6/4. 

Setting v^ = Umo^o) ^^ S^t the desired result. D 

The next four lemmas are useful to pass to the thermodynamic limit in the Coulomb 
term (Lemmas 14.31 [¥^ and [¥3|) and in the kinetic energy term (Lemma 14. 6p . 



Lemma 4.3. There exists a constant C G M+ such that for all L G N*, 

VPL e L^er(rL) H L^^A^l), D^Apl, Pl) < C (llPLllii^^(r^) + IIPilli6/5(r^J 

Proof. It is well-known (see e.g. [5]) that 

Vx G Ti, Gmix) = \x\~^ +g{x), 
with g e L°°{Ti), and that for ah L G N*, 
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Let X = {i? G 7^l I 3(x, y) G Ti X Ti s.t. x-y = R}. It holds 

y{x,y)eTLxTL, < GTZr^ix - y) < J^l^ - V - LR\'^ + L-^Wgh^- 

Therefore, for all L G N*, 

DnL{PL,PL) = GnL{x-y)pL{x)pL{y)dxdy 



XTLix)\ pL{x)\xTL{y) \pL{y)\ 

X — y — LR\ 



^ ST X^rL\X)\pL\x)\KVL\y)\PL\y)\ ^j^j^, , r-l|l„ll lU l|2 

< 2^ / / \ r^i dxdy + L 5 l°° /?l Lli (Fr) 

■^^_ /1D.3 /1D.3 \X — V — LR\ -^perl-i L) 



Rex 

< C" 1 1 XFi PL 1 1^6/5(83) + \\9\\l°°\\pl\\Ii^^(t^) 

= C'WplW 6/5,„ . + \\g\\L^\\PL\\Li (Fr)' 

where C is a constant independent of L and p^. Let Ci be the Sobolev constant such that 

Vui G -??per(ri), ||wi||i6^^.(ri) < Ci||vi||^i^^(r,). 

By an elementary scaling argument, it is easy to check that the inequality 

VvL G hI^,{Tl), lbL||L6,,(ri) < Ci\\vL\\H^^^{ri^) 
holds for all L G N*. Thus, for all vl G -ffpgj,(r/,), we obtain 

which completes the proof of Lemma 14.31 D 



Lemma 4.4. Let u G L-'^(M'^) n L^(M'^) and vl G Li^^^iVL) defined by vl\tl = ^Ir^ /or a// 



LgN*. Then 



Ymi D'ji^{vl,i^l) = D{v,v). (36) 

L— >oo 



Proof. Let ^1 := |ri| / Gi and L^ be the first Brillouin zone of the lattice TZl (that is 

the Voronoi cell of the origin in the dual space). Note that 7^2 = L^^Ti-i and F^ = L^^LJ. 
Let K > 0. We have 

2 



^•^^^ ^ fceL-l7^•\{o} ' ' 



fceL-l7^*\{o} 



LI 

feeB^nL-iTe^iMo} 

ICfe.Ll^^LJ'^ 



+4vr J: ^:^:|^, (37) 



1^ 



where Bk is the ball of radius K centered at 0, B'j^ = W^ \ Bx, 

CkA^L) = \rL\-'/' [ UL{x)e-'''-^ dx, 



and 



Cfc 



A'^l) = \n\-'^'ck,L{^L) 



1 



i/(x)e ''^'"'dx. 



As u £ L'^ 



(27r)3/2 
, \ck,L{i^L)\ < {27^)-"'/'^M\L^R3) for all k and L, P G L°°{R^), and 

Clearly the first term in the right hand side of (j37p goes to zero when L goes to infinity. 
Besides, 

|p(t)p 



E 

fceBxnL-l7^*\{o} 



|p^| |cfc,L(t^L) 

I ^1 \k\^ 



L^^ -IBk \k? 



dk. 



Lastly, 



\CkMVL) 



1^ Ifcp 

A:e-B^nL-i7?.*\{0} ' ' 



< 



< 



E 



^fceB^,nL-i-R,*\{o} 
1 



|Cfc.L('^L)|^ 



1/2 



1/2 



V ir*i — 

^ ^ \fcgBj,nL-i7e5;\{o} ' ' 



X] |cfc,L(i^L) 

^A:eBf^,nL-i7e*\{0} 
1/2 

ll^llLifM3)||'^||L2 



1 



L^oo (2^2^^)1/2 ' 

It is then easy to conclude that (|36p holds true. 



^\\L^m?)\\^\\L^ 



D 



Lemma 4.5. Let {pl)l£W be a sequence of functions of L'^^^{]R.^) such that 

1. for each LeW, pl£ -^pcr(rL); 

2. there exists C G E+ such that for all L G N*, 

PL <C and Dtil{pl,Pl) < C; 

3. there exists p G D'(M'^) such that {pl)l€N* converges to p in ^'(M^). 
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Then p & C and 



D{p,p) <limmf DnApL,PL)- (38) 

L—^oo 



In addition, for any p > 6/5 and any sequence (fL)LgN* of functions o/L|'^^(M^) such that 
vl G -^per(r_L) for all L G N*, which weakly converges to some v G L^q^(M^) in L^^^{M?), it 
holds 

V0 G Cr(M3), lim Dn^{pL,VL<P) = D{pM)- (39) 



Proof. Let Wl the unique solution in H'^^^iTL) to 



-AWl = ^t^(pl-\Tl\-^ I pi\ 



It holds 



Wl TeL-periodic, / Wl = Q. 



i-^ \VWl\^ = Dn,ipL,PL) - 9iL-' (^1^ Pl) < C, 



(40) 



(41) 



where gi := \Ti\ / G-ji^ > 0. Hence the sequence (||VM^l||j;^2 (ri^))L€N* is bounded. 
Jti "" 

By Sobolev and Poincare-Wirtinger inequalities, we have 

and by a scaling argument, we obtain that for all L G N*, 

VFi G Fi,r(rL) s.t. /" yz. = 0, ||Vi||i6^^(r,)<C;||VVi|U2^^(r,), 

where the constant C[ does not depend on L. Thus, the sequence (||W^l||l6 (rj^))LeN* is 
bounded. Let C G M+ such that 

VLGN*, ||H^L||L«,,(r,)<C and ||VH^i|U2^^(p,) < C, 

and let {Rn)n&n be an increasing sequence of positive real numbers such that lim„_j.oo Rn = 
oo. Let i? > 0. For L > 2R, 

||H^l||l6(b^) < \\WL\\Ll^^iT^) < C and \\VWL\\LHBn) < W^WLhl^^iT,) < C. 
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We can therefore extract from (VFL)LgN* a subsequence {W]^o)neN such that {Wio {bh )neN 
converges weakly in H^(Bjig), strongly in LP{Bjig) for all 1 < p < 6, and almost everywhere 
in Bfifj to some W^ G H^(Bji^^), for which 

\\W'\\LHBno)^C and \\VW^h2^B^j<C. 

By recursion, we then extract from (Wik )„gN a subsequence (W^k+i)ne'N such that (W^k+i \bji )neN 
converges weakly in H^[Bfi^^^), strongly in LP{Bji^^-^) for all 1 < p < 6, and almost every- 
where in Bji^^-^ to some W^^ G H^{Bji^^-^), for which 



rk+l\ 



rk+l\ 



\W^'"\\LHBn,,,)<C and ||Vt^«+l^.(^,^^^) < C. 



^A:+l 



LJ^ and let W be the function of H^^^ 



Necessarily, W+'Ibi,^ = W\ Let U 

by M^lsjj = VF'^ for all A: G N (this definition is consistent since W^~^^\b[, 



(42) 

defined 

= W^). 
for all 



The sequence (WL^)„gN converges to W weakly in H^^^i^), strongly in L^^^ 

1 < p < 6 and almost everywhere in R^. Besides, as (j42p holds for all A;, we also have 

||W^||l6(r3) < C and ||Viy||i2(K3) < C. 

Letting n go to infinity in (|40p with L = Ln, we get 

-AW = A-Kp. 

Introducing the dual 

C' = {V e L^{M.^) I VV G (l2(r3))3} , 

of C, we can reformulate the above results as VF G C and — Aiy = Anp. As —A is an 

isomorphism from C' to C, we necessarily have p & C. From (|41l) . we infer that for each 

R>0, 

1 



47r 



|VW"||i2(B^) < liminf D7ei(pL,/>L 



L— >oo 



Letting i? go to infinity, we end up with (|38|) . By uniqueness of the limit, the whole sequence 
(W^L)LgN* converges to W weakly in H^^^{]^^), and strongly in Lj^^(M'^) for all 1 < p < 6. 

Let p > 6/5, {vL)LeN be a sequence of functions on L^^^(M'^) such that vl G Lpc^i^L) 



for all L G N*, and converging to some v G Lif^^ 
We have, for L large enough. 



weakly in Lf 



loc 



and 



-'pcvy 

gC- 



DnLiPL,VL(, 



WiVLc/i-giL ^ 
(Wl(I))vl - giL 



Supp(0) 




VL(, 



Supp{<f>) 



L— >-oo 



Supp(0) 



W(l)v = D{p,v(l)), 



which proves (|39|) . 



D 
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Let us introduce for each L G N* the bounded hnear operator 

iL : L\^^) ^ L%,[Tl) (43) 

and its adjoint i^ G C{Lp^^{r l) , L'^ {M.^)) . Note that for all vl G L'^^j.{Tl), i\vL = Xr^^L 
and zl«2 = lL2^^(rL)- As C^(]R^) C i/^(M^), the domain of the self-adjoint operator 
(i^pgj. — Ep)^'^, the function (-ffpgr — Cp)^ '^^ ^^ ^'^ L^(M^). Using the same abuse of notation 
as above, we can also consider H^^^ as a self-adjoint operator on L'^^^^{Tl) with domain 
H'^^^{Tl) and introduce the function i*^{H^^^ — e^Y^'^iL't'^ which is well-defined in L^(R^). 

Lemma 4.6. Let G C^(M'^). The sequence {i*j^{Hp^j. — ep)-'^/^*^^)^,^^. converges to 

Proof. According to Bloch-Floquet theory |18) . each / G L'^{M.^) can be decomposed as 

fix) = -^ [ Mx)e^'-^dk 

where fk is the function of Lpgj.(ri) defined for almost all /c G M^ by 



Recall that 

y{f,g) G L2(r3) X L^R^), (/, 5)^2(^3) = -i- j Uk,9k)Ll^^iv,)dk. 

The operator -ffpgrj considered as a self-adjoint operator on L^(M'^), commutes with the 
translations of the lattice IZi and can therefore be decomposed as 



^per = TT^ / (-f^por)fc dk 



["P* I / 

mi Jti 

where {H^^^)k is the self-adjoint operator on L'i^^^Y'i) with domain H^^^iJ^i) defined by 
{HlJ, = -A - 2^fc • V + |A;p + ^(p^.J^/s + ^^0^^, 



Let i;^ and f/' be two functions of C^(M ). Simple calculations show that for L large enough 

A;er*n-R,£ 
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and 

||il(ff0,, - 60)V2^^<^||2^ = \\(H' - ,0)1/2^||2^ (45) 



per ^tJ "^y^\\L^(\K-i) ll\""per 

The sequence (i2(-ffper ~ (-p)^ iL4')LeN* therefore is bounded in L^(R^), hence possesses a 
weakly converging subsequence. 

Besides, the function k i— )• ((-f^per ~ ^f)a; 4'k^'4'k)L'^ (Vi) is continuous on T\ since 

((-f^per-4)fc (t>k,'>Pk)Ll,^{Ti) = ((^per-eF + l)fc^(-f^per-eF)fc (^fc, (-f^per -4 + l)fc^fc)L2^,(ri) 

with k ^^ (f)k and k i— )• {H^^^ ~ ^f + ^)k4'k continuous from T\ to L'^^^iTi) and k i— )• 
(if°^, - e^ + l)-^(ijOg^ - e^)^''^ continuous from T\ to £(L2^j,(ri)). Interpreting (glD as a 
Riemann sum, we obtain 

The above result allows to identify (^pcj. — ep)^'^(/> as the weak limit of the sequence 
(z2(-ffper ~ ^y)^ '^L'P)Lefi*i a-^^ (|^5|) shows that the convergence actually holds strongly in 



4.2 Proof of Proposition 12.11 

Let (i'n)neN be a minimizing sequence for ([2]). As each of the three terms of E'^^ {p^^'^ , ■) 
is non-negative, the sequence (un)neN is clearly bounded in ffpgj.(r), hence converges, up to 
extraction, to some u S i^pg^ (F), weakly in ifpgj.(r), strongly in Lper(r) for each 1 < p < 6 
and almost everywhere in R^. Passing to the liminf in the energy and to the limit in 
the constraint, we obtain that u satisfies E'^^ {p^^^,u) < I'ji{p^^^,Q) and /p u^ = Q. 
Therefore, -u is a minimizer of Q. As \u\ G Hl^^{T), ^^FW^^nuc^ |^|) ^ E^TFW^^nuc^^) ^^^ 

/p |up = /pU^, |n| also is a minimizer of ([2]). Up to replacing u with |u|, we can therefore 
assume that n > in M^. Clearly, —u also is a minimizer of ([2]). 

Working on the Euler equation ([3]), we obtain by elementary elliptic regularity argu- 
ments [9] that u € -f^per(r) "-> C^(M^) n L°°(M3), and it fohows from Harnack's inequality 
|9] that u > in R^. 

Lastly, vq is a minimizer of ([2]) if and only if po = f q is a minimizer of 

mi{Sl^'^ip^^'^,p),pGJCn,Q}, (46) 

where 

4^W(^nuc^ ^) = Cw ^ I VVpI' + Ctf ^ p'/' + ^Vnip^''' - p, p^""' - p), 
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and 

iCn,Q = |p > 0, Vp G ^pcr(r), J^p = q}- 

The functional p i— )• £^ {p'^^'^,p) being strictly convex on the convex set /C, (I46p has 
a unique minimizer po and it holds p^ = v? > 0. Any minimizer vq of ([2]) satisfying 
Vq = Po > 0, the only mininiizers of ([2]) are u and —u. 

4.3 Existence of a minimizer to (1151) 

The existence of a minimizer to (jlSp is an obvious consequence of the following lemma. 



Lemma 4.7. It holds 

3/3 >0 s.t. Vz^gC, V«gQ+, /3|b||Hi(R3)<^"(«), (47) 



Vi/GC, V«gQ+, \\ul,M\c<£"{v) + \\vYc + \Mh (48) 

and /or eac/i u ^ C, the functional E^ is weakly lower semicontiuous in the closed convex 
subset Q+ of Q. 

Indeed, if (fn)neN is a minimizing sequence for (jlSp . we infer from (j47p and (|48p that 
(f„)„gN is bounded in Q. We can therefore extract from (un)neN a subsequence {vnf,)k&f>i 
weakly converging in Q to some v^ £ Q. As Q_|_ is convex and strongly closed in Q, it is 
weakly closed in Q. Hence v^ G Q+. Besides, if being weakly l.s.c. in Q_|_, we obtain 

.S^(t;^) <liminf£:'^(?;„J = r. 

fc— >oo 

Therefore Ui, is a minimizer of (fTSl). 



Proof of Lemma \4-7\ Using (jQ]), pT]) . Lemma |4.H and the non-negativity of D, we obtain 
that for all i^ £ C and all v G Q+, 



f'^(„) > ^m^/^\\v\\l2 



2 
3 

and 



£'{v) > l|V.||i.(^3)-gMV^ + ||Fp° 
Therefore, there exists some constant /3 > such that 
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Besides, for all z^ G C and all v E Q+, 

1 1 

Hence (|I8|) . 

Let 7; G Q-(_ and (f„)„gN be a sequence of elements of Q+ weakly converging to v 
in Q. As (wn)neN is weakly converging, it is bounded in Q, which means that {vn)neN 
and {UpQ^Vn)n£N are bounded in //^(M^) and C respectively. We also notice that {v'^)n£N is 
bounded in L^(M.^) D L^{R^) ^ L^/^{M.^) ^ C. 

Therefore, we can extract from (fn)neN a subsequence (fnfc)fcGN such that 

• {£''{vnj)km converges to / = liminf„_^oo'?'^(t'n) in M+; 

• {vn^.)k&N converges to some v G H^(W^) weakly in /7^(R^), strongly in L^^^(IR^) for 
all 1 < p < 6 and almost everywhere in M'^; 

• (^pcr^nfc)fceN weakly converges in C to some w £ C; 

• (f^ )fcgN weakly converges in C to some z £ C. 
We can rewrite the last two items above as 

VFgC, / ut,Vn,V ^ I wV, and / <F ^ / zV. 

Together with the strong convergence of {vnf,)keN to v in L^q^(]R^), this leads to u^^^v = w £ 
C and z = v^. This in turn implies that {vnf.)keN weakly converges in Q to v. Therefore 
V = V. Finally, {vn,.)k<=N converges to v weakly in H^{M.^) and almost everywhere in M^ and 
{2Upg^Vnf, + v^f, — i^)keN weakly converges to 2Upgj.v + v^ — v vnC. 
It follows from dill) that 



((-f^per - eF)f^>^)H-i(M3),Hl(IR3) < lim inf ((Fp<,j, - eY)Vn^,Vn^) H-i{M.3),m 

By Fatou's Lemma, 

< hminf 1^^ (\ul, + .„JW3 _ |^o^^|io/3 _ ^|«o^//3(2^o^^,„^ + y^)^ . 

Lastly, 
D{2ul^^v + v'^ -V, 2ul^^v + t>2 - I/) < lim inf Z?(2M°^^t;„^ + vl^ - v, 2ul^^Vn^ + vl^ - v). 



fc— >-oo 
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Consequently, 



E^iv) < liminf £:^(v„, ) = liminf £:'^(w„ 

A;-s>oo n—>-co 



which proves that S'^ is weakly l.s.c. in Q^. D 

Clearly, the functional S'^ is C^ in Q and it holds 



V/iGQ, {r'{v),h)Q>,Q = 2[{{H^,,-e^p)v,h)H^i^^3)^H^ 



5 



+ 3 Lil^P- + 



,0 



y\V3 _ I |7/3 _ I |4/3 \ ^ 
u\ I "perl I "perl '^ J "' 



+D{2u%,v + V^ -V, (n°er + v)h)\. 



The minimization set Q+ being convex, v^, satisfies the Euler equation 

Vi;GQ+, {£'"{v^),{v-v^))q,^q>{). (49) 

Let n,y = Uppj. + Ujy and 

y = ^pcr - 4 + ^1^.1'/' + (2<r^z. + Z;2 - z.) * I . |-1. 

The function Uy satisfies u^ G //^^^.(M^), n,y > in M'^, and 

= 2^£'"iv^),{Vu + (t>-Vu))Q',Q- 

This imphes that for all £ C;?°(R3) such that (/> > in M^ 

Vuu • V(/> + / Vuu(t> > 0, 



since Vi, + (j) £ Q+. Therefore, Ui, is a non-negative supersolution of —Au + Vu = 0, with 
y G L[^^ (M^). It follows from Harnack's inequality (see Theorem 5.2 of |20| ) that either 
u,^ is identically equal to zero in R^, or for each bounded domain Q of M^, there exists r/ > 
such that Vu > — ^per + ?/ in f]. As the first case is excluded since — wSer ^ 2+) (|49p implies 
£^'{vp) = 0, which means that Vj/ is a solution in Q+ to the elliptic equation (jlSp . 



Remarking that 



f^(t;.)<^^(0) = ^z?(z.,^) = i||z.||2, 



and using (|47p . (|48p and Lemma [4.31 we finally get the estimate (|17p . 
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4.4 Uniqueness of the minimizer to ( llSp 

Noticing that 

Q+ = {ve H\R') I {u%, + vf - p%, G C, u%, +v>0}, 
we obtain that u* is a minimizer to (|15p if and only if p^, = (Upgj. + w*)'^ is a minimizer to 

inf{g(/)), pe/C} (50) 

where 

Gip) = J{p) + j^^ [p'/' - {pW" - \{pWHp - pU) + \d{p - p\ 

and 

To see that /C is convex and that Q is strictly convex on /C, we first introduce the set 



per - 2^, P - Pper " ^) 



where £"'(1^ ) denotes the space of the compactly supported distributions, and observe that 
for all p & IC, 

J{p) = j^^ (iVVpP - iVnO,/ + {^-{pl^fl' + Fp^er - 4) {P - P' 



9° ^ 
^per/ 

Reasoning as in the proof of the convexity of the functional p i— )• L3 | V-^/pp on the convex 
set {p > I yfp € i?^(IR^)} (see e.g. |15|). we obtain that K, is convex and that J is convex 

on /C. It then follows that Q is strictly convex on /C. We finally conclude by a density 
argument. 

As Q is strictly convex on the convex set /C, (150p has at most one solution. Therefore, 
Py = (Upgj, + Vi,)'^ is the unique solution to (j50p . and fj, is the unique solution to (jlSp . 

4.5 Properties of the unique minimizer of (I15[) 

The Euler equation (jlSp can be rewritten as 

- ^V, + K<r = / + (Z^ * I • r'Xr, (51) 

where 

/ = (4 - ^per)^. - \ (kper + ^^.l'^' " WW'^) " ((2<r^. + Z;^ _ i.) ^ | . |-1) ^^, 
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and where V^ = (2Wpgj.fi, + ^^) * I " I ^ satisfies 

-AV, = i7r{2u%^v, + vl). (52) 

We know that v^ G H^{R^) ^ C and that K G C since (2u°gri;i, + vl) £ C. Adding up 
()5ip and (j52p , we obtain that Wi, = v^ + Vy is a solution in C' to 

- Aiy, + uOgjTy, = 7+ (i. * I . |-^)n°er, (53) 

where f = f + (87r+ l)v^^j.Vu + ^i^v^ G L^(R'^). Since Up^j. satisfies Q, the ehiptic equation 

-Aw + Uper^ = / 



per 

2 



has a unique variational solution in H^{M?), which we denote by Wy. Clearly Wy (L H 
The function Wy = Wy — Wu £ C then is solution to 

- AWu + Up^^Wu = (l/ • I • r^)liper- (^4) 

Introducing p^, = — (47r)~^A-u;j^ € C, (I54p also reads 



Therefore, 



47r-^ = {u- py) -k 

^per 



47r / -^ = D{u -py,p„) <oo, 

Jr3 «^er 

which proves that /)j, G L^(]R^), hence that {ly -pu)*\- \~^ G L^(M^). As 

(z^-(2Upgj,v,,+t;^))*|-|"^ = i/*|-|"^-K = (z^-/?zy)*|-r"^+'u;^-K = (z^-/>,y)*|-|"-^+Vi.-i(;i., 

we obtain 

1>l = {u- {2ul,,v, + vD) * I • r^ G L^W"). 

Introducing p^ = u — {2Upgj.v^ + v^), the above statement reads 

Jr3 |fc|4 

Therefore, 

, .. , , , \l/2 / /. ir0^l.M2 

-Or 



^l^j^o,,,,,, , _L(/^j,...)"7|^^^.. 



1/2 



12 \ 1/2 



Lastly, rewritting (|5ip as 
we conclude that Vi, G -fT^ 






-Av, = f + <^lu%„ 



28 



4.6 End of the proof of Theorem 13.11 

We have proven in the previous two sections that: 
1. (jlSp has a unique minimizer v,^; 



2. if {vn)n£'N is a minimizing sequence for (jlSp . we can extract from (t^n)nGN a subse- 
quence {vn^)keN which converges to v^,, weakly in H^(R.^), and strongly in L[^^(M'^' 
for all 1 < p < 6, and such that {u^^j.Vn^.) keN converges to Up^j.v,y weakly in C. 

By uniqueness of the limit, this implies that any minimizing sequence (vn)nGN for (fT5] 



converges to v^,, weakly in H^(M^^), and strongly in L|'q^(M'^) for all 1 < p < 6, and that 
{Up^j.Vn)neN Converges weakly to u^^^Vp in C. Lastly, the existence of a minimizing sequence 
for (jlSp satisfying (|20|) is a straightforward consequence of Lemma 14.21 

4.7 Thermodynamic limit with a charge constraint 

Let u G L"'^(IR^)nL^(M^). Clearly, Vi/,g,L is a minimizer to (I27p if and only if tipgj,+Ui,^g^i is a 
minimizer to ([2]) with 7^ = TZi, p™'^ = Pp^^ + i^L and Q = ZLi^ + q such that liper + ^i^.g.L ^ 
in R^. It follows from Proposition 12. II that (j27p has a unique minimizer Vy^q^L, which satisfies 
Vu,q,L e H^cA^l) ^ C2(M3) n L°°(IR3) and u^^^ + v^^q^i > in M^, and the Euler equation 
(f29|) for some ^^.^^^l G M. 

Let a = iril^"*^ Jp "Uper- For L large enough, a^ + g/lFil > and the constant function 
zl = —a + y^^+^/ITlI satisfies z^ > — ^per everywhere in M.^ and 

/ {2u%,ZL + zl) = q. 

Using Lemma |4.H Lemma 14.41 and the fact that I-Zl] < CL^^ for some constant C inden- 
pendent of L, we obtain 

£1[vu,,,l) < £1{zl) 



,,0 I zr l^°/3 _ L.O |10/3 _ 12u,0 |7/3 
"per T^ ^-L I I "per I q I ""per I 



+ 2-^^L (2^iper^L + zl- ^L, ^u^^.zl + zl - Ul) £^ D{u, v). (55) 
Besides, reasoning as in Section [4.31 we obtain 

where the constant /3 > is the same as in (|47p . and 

1 
2' 

< 8l{vL) + DnMA) + DnA^L,yL)- (57) 



^VL^Q+^L, Dti^{u1^^vl,u°^^vl) < £1{vl) + -Dti^{vI - VL,vl - vl) 
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We infer from ([55|) and (|56p that for each g G M, there exists Cq £ M-|- such that 

VLGN*, \\vu,,,L\\H^^jr,)<Cg. (58) 

By a diagonal extraction process similar to the one used in the proof of Lemma [4. 51 we can 
extract from (fv,(j,L)LgN* a subsequence (wj/,ij,Lj.)fcGN which converges to some u^, G if^(R'^), 
weakly in Hy^^{M.^), strongly in L^^^(M^) for all 1 < p < 6 and almost everywhere in R^ and 



such that 



lim £L,ivu,q,Lk) =lircimf£Livu,q,L)- 



fc— >-CXD L— >-oo 

In particular Ui, > — Uper almost everywhere in 



Let us now prove that Up^^^u,^ G C. First, we notice that it follows from (|55p . (j57p and 



"per '^ 

Lemma 14.31 that there exists a constant Cq such that 



DnA'^lerVu,q,L,ul^rVu^q^L) < Cq 



Besides, 



'^pcr'^'^,g,L 



*per 
1 

2 



(59) 



V,g,L 



<^(M + ^.^) 



'2 



hence in the distributional 



and iUpevVi',q,Lk)keN converges to Up^^Uu strongly in L 

sense. It therefore follows from Lemma 14.51 that Up^j.Ui, G C. Thus, Ui, G Q+. 

As (j29p holds in /7~].(ri), we can take u^^^ as a test function. Using (jlOp . we obtain 



,0 



.0 



^^v,q,L\ZL + / fv,<j,L''iper I " / '^ ( I'^^per + ^i^.g.il' ' " ~ l^perl ' ' " ~ l^perT' "^!^,g,i j ^per 



.0 



|7/3 _ 1^0 |7/3 



,0 |4/3, 



Using (I58p . (I59p and Lemma [4. 11 we obtain 



/JO 



''^per + '^>^,q,L 



7/3 



,0 1 7/3 



,0 |4/3 



"perl 



"perl 



^^i.,q,L U 



,0 
"per 



< C:L3/2, 



\DTlr. {{2u%,V,,q,L + V^q^L - ^l), «r + ^..g.^Xr) | < Cql""/^ , 



Vu,q,LU 



per 



<^(M + c,^)> 



for some constant Cq independent of L, which allows us to conclude that {lJ^u,q,L)Lm* goes 
to zero when L goes to infinity. 

Note that using Lemma 14.51 we can pass to the limit in the Euler equation (|29p in the 
distributional sense, and prove that u^ satisfies 



(ttO 
v-"per 



Ef)^!^ + o ( l^iper + Uuf^^ - \U 



|7/3 _ |^0^j4/3^^ 



perl 



"perl 



+ ((2^iper'"i' + ul 



l^)*|-|"')(<r + ^.)=0. 



(60) 



30 



We are now going to prove that £^{uu) < £^{vu), which implies that Ui, = v^ and, 
by uniqueness of the Umit, that the whole sequence (fi/,g,L)LeN* converges to Vu weakly in 
Hl^^{R^), and strongly in ^[^^(M^) for all 1 < p < 6. 

Let e > 0. From Lemma [4.21 there exists v^ £ Q+ Pi C^(IR'^) such that 



and 

£''M<£''{vt^^)<£'^{v,,,) + e. 

For L large enough, the T^^-periodic function v^, ^ defined by v^ ^\r^ = vlJrj^ is in the 
minimization set of (f27|) . Using Lemma [4.41 and the fact that v^g is compactly supported, 
we have for L large enough t^^ ^ € Q+,l and 

£LiVu,g,L) < Sli'^t^q^L) = ((-f^pcr -eF)^^,<7)'^^,g)H-i(IR3),Hi(IR3) 



. / f\uO <y^ |10/3_|^0 |10/3_5|0 |4/3/2^0 
/ \ I'^per ^ '^u,q\ I "per I o I "per I \-^"'p 

Jr3 \ o 



t \2' 



1 

2' 



+ i;DTir^ (2^iper<g,L + «g,L)^ " J^L, 2n°er<g,L + «g,L)^ " ^^l) 



Therefore, for each e > 0, 
for L large enough, so that 



limsui>£l{vu,q,L) <£''{vu)- (61) 

L—^oo 



We are now going to prove that 

£''{u,)<limmf£Uvu,q,L). (62) 

For each A; G N, we denote by 

Vk ■■= il^Vu,q,Lk and Wk := ^IJ^per - 4y^'^v,,,q,L^, 
where the operator ii,, is defined by (j43p . As H'^fcllizms) = ||f;y,g,Lj. 11/^2 (p^ ) and 

||u'fc||i2(K3) = ((i?por - ^F)Vu,q,Lk,V>^,q,Lk)H-,\iTL^),H-^ArLk)' 

we can extract from ('yfc)fcGN and (wfc)fcgN subsequences {vk„)n&f>i and {wk^)nef>i which weakly 
converge in L^(R^) to some v G L^(M^) and w G L^(M^) respectively, and such that 



lim £" {viy^q^LkJ = liminf.S''(7;,,,q,L). 
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As (Wi/,q,Lj.)fcgN converges to u^, strongly in L^j^^(M^), we have v = Uy. Let us now prove 
that w = (-ffper ~ ^y)^^'^'^^- -^°^ ea,c\\ (p G C^(M^), we infer from Lemma 1^6) that 

jim (il^^^ (i^o,, - e'F)'/^^L,Jk„A)L^ 






nm(?y,„,iL K,,-eO)V2, 0)^, 



= (Wy, (^per - ^f) 0)l2(IR3) = ((//per " Cp) ""v, 0)l2(]r3). 

As a consequence, u; = (-f/^per — Cp)^ ^^j^- 

Using the weak convergence of Wk„ to w = (^per ~ ^fY ^^^ Fatou's Lemma and 
Lemma 14.51 we thus obtain 



s'-M = iKi/o -4)V2^^||2 



+ / (|n0,, + n.|i°/3 



^perl ol'^^perl [^'^per'^^i' ~^ ''^uJ j 



< hminf£:''(7;,^„L ) = hminf .^''(u^^l). 

n->oo " L— i-oo 

Hence (|62p . Gathering (j6ip and (j62p . we obtain that £^{ui,) < £^{vy) and therefore that 
^i/ = Vy since Ui, G Q+ and (llSp has a unique minimizer. 



4.8 Thermodynamic limit without a charge constraint 

Let {vn)n&i be a minimizing sequence for (|28p . For all ?? > 0, for n large enough, 

Thus, (un)neN is bounded in H^^^{T]^). Extracting a converging subsequence and passing 
to the liminf in the energy, we obtain a solution Vi^^l to (128]) . such that 

/5|l^^,ill^i„(ri) < ^DTi^i^L^t'L)- (63) 

We also get 

DnL{ulerVu,L,ul^,Vt,^L) < C, (64) 

for some constant C independent of L. 

Clearly, u^^^ + Vi^^l is a non-negative solution to 

inf {Ell'^ipl-^ + J.l,Wl), Wl G ^^er(rL)} • 
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Reasoning as in the proof of Proposition 12.11 we obtain that Upgj. + Vy^L is the only non- 
negative solution to the above problem, and therefore that v^^l is the unique solution 
to (|28p . Besides, v^^l G -f^perl^L), "Uper + '^u,L > in R^, and v^^l is solution to the Euler 
equation (j30p . which holds in //"{.(r^^). Taking u^^^ as a test function, we get 

/I I |7/3 I |7/3 I |4/3 \ 

o M'^^per ~r Vu^l\ ~ l^perl ~ Pperl ^1^,^) ^per 

+^7ei ((2n°erf^i',i + '^l,L " ^^l), f^i^,L^iper) + ^Te^ {{2ul^^Vu,L + W^,L " Z^l)> (^per)^) = 0- 

We now remark that the third term can be rewritten as 

+ / {2U%V,,L + <L - i^iXer, (65) 

where, as above, gi = |ri|^^ Jp G-ji^ and where VFpgj. is the unique solution in H'^^^iVi) to 



VF^g, 7^l-periodic, / TV V = 0- 
'ri 



We finally obtain 

As the right hand side is bounded by CL^'"^ for a constant C independent of L, it holds 

lim {v- {2ul^^Vy^L + vIl)) = 0. 



Proceeding mutatis mutandis as in the previous section, it can be shown that the se- 
quence (ujy,L)LGN* converges weakly in /7[|^^(]R^) and strongly in L]'^^(M'^) for all 1 < p < 6, 
towards the unique solution Vy to (jlSp . 
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